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Abstract. In the reduced phase space of electromagnetism, the generator of 
duality rotations in the usual Poisson bracket is shown to generate Maxwell's 
equations in a second, much simpler Poisson bracket. This gives rise to a hierar- 
chy of bi-Hamiltonian evolution equations in the standard way. The result can be 
extended to linearized Yang-Mills theory, linearized gravity and massless higher 
spin gauge fields. 
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1 Introduction 

A cornerstone of soliton theory is the discovery that the evolution equations are Hamilto- 
nian systems 03121. hi this context, the occurrence of hierarchies of evolution equations 
sharing the same infinite set of conservation laws can be understood as a consequence 
of the existence of a second, compatible Hamiltonian structure giving rise to the same 
evolution equations (31 0. 

The equations of motion associated to the theories for the known fundamental forces 
of nature, electromagnetism, Yang-Mills theories and gravitation, are variational and thus 
Hamiltonian. This is no coincidence, since these theories are fundamentally quantum, at 
least the first three of them, and only for variational theories quantization is sufficiently 
well understood. 

In order to have Poincare invariance, respectively diffeomorphism invariance, mani- 
festly realized, most modern investigations of these equations are carried out in the La- 
grangian framework. This could be the reason why the bi-Hamiltonian structure underly- 
ing these equations and discussed below has hitherto remained unnoticed. 

At the heart of our analysis is an important exception to this paradigm, namely the 
question whether the duality invariance of the four dimensional Maxwell or linearized 
gravity equations admits a canonical generator. This question has been answered to the 
affirmative in the reduced phase space of these theories and generalized to massless higher 
spin gauge fields [El [6113. 

We will show in this letter that the reduced phase space formulation of massless higher 
spin gauge fields is bi-Hamiltonian. The second Poisson bracket on reduced phase space 
turns out to be more natural than the one induced from the covariant variational principle, 
while the generator for duality rotations plays the role of the second Hamiltonian. 

We will start with the most familiar case of Maxwell's equations, for which all details 
will be given. This result trivially generalizes to Yang-Mills theory with an invariant, 
non degenerate metric, linearized around a zero potential by decorating the expressions 
obtained in the electromagnetic case with an additional Lie algebra index. 

The reduced phase space of linearized gravity around flat spacetime consists of two 
symmetric transverse and traceless potentials [H. The entire analysis of the electromag- 
netic case then carries over in a straightforward way by taking care of the additional 
spatial index. The same goes for massless gauge fields of spin higher than 2. 

Several generalizations and extensions are suggested by the result reported here. A 
first exercise consists in studying the consequences for symmetries and conservation laws 
of both the Maxwell and the higher spin equations and compare them to known results 
(see e.g. [[H [10l [CD and references therein). Another obvious question is to investigate 
more general backgrounds. For instance, the generalization to massless spins propa- 
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gating on (anti-) de Sitter spaces instead of Minkowski spacetime and the inclusion of 
fermionic gauge fields should be straightforward. In Yang-Mills theories (anti) self-dual 
backgrounds could be promising in view of their close connection to integrable systems. 

The most important problem is however the inclusion of interactions. When compar- 
ing to the Korteweg-de Vries equation for instance, the present work corresponds to the 
bi-Hamiltonian structure for the linearized equation. The question is then to find inter- 
actions that preserve this structure. This will most probably not work for the standard 
interactions, as it is known that duality invariance does not survive the standard Yang- 
Mills type coupling, nor the extension from linearized to full gravity (51 [13. A promising 
candidate for an interacting model in which duality invariance does survive is Born-Infeld 
theory^. We plan to address at least some of these questions elsewhere. 

2 Electromagnetism 

The first order Hamiltonian action for free electromagnetism is 



Here, Aj is the vector potential, E % the electric field and A the Lagrange multiplier for the 
Gauss constraint. We use x° = t and O is the curl, (OA) 1 = e ljk djA k . Unless otherwise 
specified, the summation convention is used, Roman indices i,j,k... take values from 
1 to 3, are raised and lowered with the Kronecker delta, while is completely skew- 
symmetric with ei23 = I. 

The decomposition of A4, E % into transverse and longitudinal vectors Ai = Af + 
Af, E l = E T% + E Ll allows one to eliminate the pure gauge degrees of freedom Af , E lL 
associated with the Gauss constraint diE % = 0. By introducing a second vector potential 
Zi such that E lT = (OZ) % the reduced phase space action has been shown to be invariant 
under duality rotations 0. This is most transparent in the double potential notation lfi"3l 
where A" = (Ai, Zi). Roman indices a,b,c . . . take values from 1 to 2, are raised and 
lowered with the Kronecker delta, while e ab is completely skew- symmetric with e 12 = 1. 
In this notation, the reduced phase space action can be written as 




(2.2) 



(2.1) 




(2.3) 



(2.4) 



'The authors thank X. Bekaert for pointing this out. 
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where in the second expression for the Hamiltonian, we have used that O is "self-adjoint", 

J d 3 x gi {OfY = j d 3 x {Og)d\ (2.5) 

and 

(0 2 f T Y = -Af Ti , (2.6) 

with A = did 1 the Laplacian in fiat space. The standard Poisson bracket determined by 
the kinetic term is 

{Aj a (x) J A Tb3 (y)} 1 = e o6 A-V^J (3) (a; - y) = e" 6 A" 1 (O y 8f (3) (x - y)) j , (2.7) 

where 8jj(x — y) is the transverse delta function, see e.g. Ifl4l section Ai.2. In vacuum, 
Maxwell's equations for the physical degrees of freedom read 

d A Tai (x) = {A Tai (x),H 1 } 1 = -e ab (OAlY(x), (2.8) 

while the generator for duality rotations is 

H ° = ~IJ d3x A T(OA T a ) u {H , = 0. (2.9) 

Remarkably, the duality generator is simply an SO (2) Chern-Simons action [Tl5l . 

When presented in this way, the second Hamiltonian structure is obvious and a lot 
simpler than the one induced from the covariant action principle. Indeed, a natural Pois- 
son bracket on reduced phase space is simply 

{Af a (x), Aj b (y)}o = e ah ^\x - y), (2.10) 
in terms of which the duality generator is the new Hamiltonian for Maxwell's equations, 

{A^ix^Hxh = {A Tai (x),H } . (2.11) 

This is our main result. 

At this stage, one can pause and ask whether electromagnetism and its quantization 
should not be based on this new Hamiltonian structure. A good reason to favour the old, 
more complicated structure is that, by construction, the Poincare and conformal symme- 
tries admit canonical generators for the old structure, while not all of them do for the new 
one. We plan to return to this question in detail elsewhere. 

The reason why we are interested in the new structure is not so much in order to 
use it as an alternative for quantization, but because of what it tells us about Maxwell's 
equations themselves. For instance, using the bi-Hamiltonian structure, there is a standard 
way to generate an infinite set of generators that commute with the Hamiltonian and are 
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in involution, which we now briefly describe (see e.g. QT6J], chapter 7.3 for an elementary 
introduction). 

The recursion operator is obtained by contracting the new Poisson structure with the 
inverse of the original one, 

K ai bj = -5 a b (0)). (2.12) 
Consider, for p ^ 1, Kj ai = 'e al '(O* 'Ajf , or equivalently, 

Kl a Ux) = {-) n+1 e ab A n {OA T b )\x), Kl a l 2 (x) = (-f+V^+Uf (x), (2.13) 

for n ^ 0. The evolution equations of the hierarchy 

d A Tai (x) = K^ ai (x), Vp^l, (2.14) 

are also bi-Hamiltonian, 

K™(x) = {A Tai (x), Hp]i = {A Tai (x), H p -\}q, (2.15) 
where H p ^ = ( -f j d 3 x Aj a (O p A T )i, 

(_)n+l r (_\n+l r 

H 2n = J d 3 xA Tai A n {OA T a ) h H 2n+1 = J d 3 xAj a A n+1 A T ai . (2.16) 

with Hamiltonians that are in involution, 

{H n ,H m } 1 = = {H n ,H m } , Wn.m^O. (2.17) 



3 Linearized gravity 

The Hamiltonian formulation of general relativity linearized around flat spacetime is 
based on the first order action principle (SI 



(3.1) 



Sp F [h m n, n mn , n m , n} = Jdt[J d 3 x {ir mn h mn - n m H m - nU) - H PF 

H PF [h mn ,n mn ] = J d 3 x(n mn n mn - ±tt 2 + ±d r h mn d r h rm - 

- l -d m h mn d r h rn + \d m hd n h mn - \d m hd m h) , (3.2) 
H m = -2d n 7T mn , H ± = Ah - d m d n h mn . (3.3) 

Here, h = h m m , n = n m m and the linearized 4 metric is reconstructed using /i o = — 2n 
and h 0i = rii. 

Duality invariance of massless spin 2 fields has been uncovered by introducing suit- 
able potentials for the linearized three metric and its momentum (5). We will use below 
properties of these potentials discussed in more details in lfT7l l7l[T8ll. 
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Symmetric rank two tensors 4> mn decompose as 



/T + 0L + ^n, (3-4) 



'ran ^rmn 1 rmn 1 rmn) 

ran 



dm^n + d n 1p m , (3.5) 

f>L = ^(5 mn A-a m a n )^ T . 0.6) 



The tensor 4> mn contains the trace of the transverse part of <f) mn and only one independent 
component. Three independent components are encoded in the longitudinal part <^ n , 
while the transverse-traceless part (fZ containing two independent components is defined 
as the remainder, 

^mn 4>mn 0mn < Pmn.- 

(3.7) 

The elements of the decomposition are orthogonal under integration if boundary terms 
can be neglected, 

dX(p (f mn = dx{(p fmn + ( P ( Pmn + ( P <Pmn) ■ ( 3 - 8 ) 

The generalized curl |j7l [T2J is defined through, 
It satisfies 

(OWL-C (3-10) 

and is self-adjoint, 

J d*x {Of T ) mn V T J n = J d*xf T ™ iO^) mn . (3.11) 
The reduced phase space action for linearized gravity is 



S H = dt 



d 3 X 7T TTmn doh T J n - H X 



(3.12) 



Hi = j d 3 x (n^nZ + \d T h T J n d r h™) . (3.13) 

In the reduced phase space, one can make the change of variables, 

h T J n = 2 (OH 1TT ) mn , nZ = -&H 2 J n T - (3.14) 
so that the reduced phase space action becomes 

S R [H T Jn] = jdt[jd?x e ab A {OH TTa ) mn d HZ b ~ Hi , 



(3.15) 



H t = d?x H TTamn A 2 Hjr n ). (3.16) 
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The standard Poisson bracket determined by the kinetic term is 

{HlZ a (x), H T ™\y) }l = ^A- 2 (OTJT'O* ~ y)f , (3.17) 

where Smn 7 "™ (x—y) denotes the projector on the transverse-traceless part of a symmetric 
rank two tensor. The duality generator is 

D = -Jd s x H TTamn A {OH T J) mn , {D, Hxh = 0. (3.18) 

The analogy with the spin 1 case can be made perfect by the change of variables, 
H^ = ^(OA—) mn , AT 1 =-^A- (3.19) 
in terms of which 

Hl = \jSx{ {OA^r (OA T a T ) mn ) =~\j d 3 x (A^AAZ n ), (3.20) 
{A m T n a {x\ A^iy)}, = ^A- 1 (CWjTXz - y)) kl , (3.21) 
H = -D = - l -J d 3 x A TTamn (OA T a T ) mn ■ (3.22) 

All formulae of section |2] below equation (12.91 ) now generalize in a straightforward 
way to massless spin 2 fields by replacing T (transverse) by TT (transverse-traceless) 
and contracting over the additional spatial index. 



4 Massless higher spin gauge fields 

The extension of these results to massless higher spin gauge fields [fP9l (see also [[20]) fol- 
lows directly from the observation that the Hamiltonian reduced phase space formulation 
of these theories merely involves additional spatial indices Q, so that all above results 
generalize in a straightforward way. 

This can be seen for instance by starting from the approach inspired from string field 
theory, where the Lagrangian action for massless higher spin gauge fields is written as the 
mean value of the BRST charge for a suitable first quantized particle model lETl l22l l23l 
(see also [|24ll25ll for further developments). In this framework the reduction of the action 
to the light-cone gauge corresponds to the elimination of BRST quartets composed of 
ghost and light-cone oscillators (see e.g. Eol l2~7ll ). In exactly the same way, the ghost, 
temporal and longitudinal oscillators form quartets that can be eliminated to yield the 
Lagrangian gauge fixed action for a massless field of spin s ^ 1, 

^[C.d = ~\j ^ ^U^o 1 ''' - (4-1) 
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where the field <pj^ 4 is real, completely symmetric, traceless and transverse, 



TT _^TT jJTTi - Q, ^ Q. ( 4 .2) 



The Hamiltonian formulation is direct, the momenta being ixJ T ■ = dnd)J T ' ■ . 

Consider then the Fock space defined by [a 1 , <r\ = 8p |0) = 0, the number operator 
A?" = a\a l , the "string field" <fi^ T {x) = ■ • • <4 S |O)0^ T is (x) and the inner product 

I xTT j,TT\ _ I / jTT „/,TT\ / j3„iTT jTTh.. 



,€ r ) = 'l 4 .r (oi'.rj'}/ = ^.r:o/ ; '..,.r . (4.3) 



With this inner product, the generalized curl 



O = ^ k a\d 3 a k (4.4) 



is again self-adjoint. Furthermore, it squares to —A inside the inner product involving 
trans verse-traceless fields, 

O 2 = 4r [ - AiV 2 + (0 ■ a))(d ■ a) + (a f ■ a f )A(a • a) + 2(9 ■ a f )iV(<9 ■ a) 
N z 

- (d • at) 2 (a • a) - ( a t • a t)(0 • a) 2 ] =► {f s T , O^J T ) = ~(£ T , ^D- (4-5) 

The change of variables making duality invariance transparent is 

€(. = ^^ = -(0^. (4-6) 
The first order reduced phase space action becomes 



S R [A TTa ] = Jdt[Jd*x \e ab {OA™)* 1 "*' d Al T X - H x 



(4.7) 



H 1 = ~ I d 3 x A™ AA'!,'" -' -. (4.8) 



2 



Again, all formulae of section |2] below equation (12.81) . including the one for the duality 
generator, suitably generalize by contracting over the additional spatial indices. 
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